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Abstract 

We study path integration on a quantum computer that performs quantum summa- 
tion. We assume that the measure of path integration is Gaussian, with the eigenvalues 
of its covariance operator of order j~ k with k > 1. For the Wiener measure occurring 
in many applications we have k = 2. We want to compute an e-approximation to path 
integrals whose integrands are at least Lipschitz. We prove: 

• Path integration on a quantum computer is tractable. 

• Path integration on a quantum computer can be solved roughly e~ l times faster 
than on a classical computer using randomization, and exponentially faster than 
on a classical computer with a worst case assurance. 

• The number of quantum queries needed to solve path integration is roughly the 
square root of the number of function values needed on a classical computer 
using randomization. More precisely, the number of quantum queries is at most 
4.22 e . Furthermore, a lower bound is obtained for the minimal number of 
quantum queries which shows that this bound cannot be significantly improved. 
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• The number of qubits is polynomial in e . Furthermore, for the Wiener measure 
the degree is 2 for Lipschitz functions, and the degree is 1 for smoother integrands. 



1 Introduction 

Although quantum computers currently exist only as prototypes in the laboratory, we believe 
it is important to study theoretical aspects of quantum computation and to investigate its 
potential power. There will be additional incentives to try to build quantum computers if it 
can be shown that there are substantial speed-ups for a variety of problems. 

To date there have been two major algorithms for discrete problems on quantum comput- 
ers that are significantly better than on classical computers: Shor's factorization and Grover's 
data search algorithms, see ||, [5], [19], [20f] . But numerous problems in science and engineering 



have continuous mathematical models. Examples include high dimensional integrals, path 
integrals, partial differential and integral equations, and continuous optimization. 

Continuous problems are usually solved numerically; they can only be solved to within 
uncertainty e. The computational complexity of these problems on classical computers is 
often known; for a recent survey see fl23| . Complexity is defined to be the minimal number 
of function values and arithmetic operations needed to solve the problem to within e. 

For many continuous problems defined on functions of d variables, the complexity in the 
worst case deterministic setting is exponential in e _1 or in d. In the latter case, the problem 
is said to suffer from the "curse of dimensionality" and is computationally intractable. For 
some continuous problems the curse of dimensionality can be vanquished by weakening the 
worst case deterministic assurance to a stochastic assurance, such as in the randomized 
setting. Monte Carlo is a prime example of an algorithm in the randomized setting. 

A start has been made toward solving continuous problems on quantum computers in 
recent papers % % |, ||, [16], [17|] . They study multivariate integration and approximation. 



The major technical tool in these papers is the quantum summation algorithm of Brassard, 
H0yer, Mosca and Tapp that is based on Grover's iterate, see || |4j]. The essence of the 
results of Heinrich and Novak, []7|, [8], [RJ, is that intractability in the worst case setting of 
multivariate integration in a Sobolev space is broken by the use of the quantum summation 
algorithm. That is, we have an exponential speed-up of quantum algorithms over deter- 
ministic algorithms with a worst case assurance. Furthermore, there is roughly a quadratic 
speed-up of quantum algorithms over randomized algorithms run on a classical computer. 

Our paper is a continuation of the idea of using quantum summation for continuous 
problems. Summation is often required for continuous problems. Algorithms such as Monte 
Carlo and Quasi- Monte Carlo are used for a variety of continuous problems and they require 
the summation of many terms. In the worst case setting, the number of terms n is often an 
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exponential function of e^ 1 . However, if we perform summation on a quantum computer, 
this is not a show stopper, since the cost of the quantum summation algorithm depends only 
logarithmically on n. Hence, as long as n is a single exponential function of e' 1 the quantum 
cost is polynomial, and the problem becomes tractable on a quantum computer. In this paper 
we show that quantum summation is a powerful tool for computing also path integrals. 

Path integrals may be viewed as integration of functions of infinitely many variables. They 
occur in many fields, including quantum physics and chemistry, differential equations, and 
financial mathematics. Efficient algorithms for approximating path integrals are therefore 
of great interest. However, and perhaps not surprisingly, path integration is intractable on a 
classical computer in the worst case setting for integrands with finite smoothness as shown in 



25fl . Fortunately, the worst case complexity of path integration is only a single exponential 
function in e~ l if the measure of path integration is Gaussian and the eigenvalues of the 
covariance operator are of order j~ k for k > 1. For the Wiener measure, which appears 
in many applications, we have k = 2. That is why when we use the quantum summation 
algorithm, path integration becomes tractable on a quantum computer. More precisely, for 
functions having smoothness r, see the precise definition of the class F r in Section 3, path 
integrals can be computed using of order 

• e^ 1 quantum queries, 

• £ _ ( fc +T( r ))/( A: ~ 1 ) loge™ 1 quantum operations, and 
. e -(i-Hy(r))/(fc-l) log £ -i qu bits. 

Here 7(1) = 1 and 7(7") = for r > 2. For the Wiener measure, we have more specific bounds, 
which we present in Theorem 2. We stress that the cost of a quantum query depends on a 
particular applications and may be very large. 

We know more precise bounds on the number of quantum queries. To explain them we 
comment on two types of errors for path integration. The first error occurs when we replace 
the original problem by finite dimensional Gaussian integration, and then the second one 
when we approximate the finite dimensional problem by a finite sum. For simplicity, we 
assume that the both errors are bounded by e/2 so that the total error is at most e. Hence, 
we need to apply the quantum summation algorithm with error e/2. To get a 5-error, the 
quantum summation algorithms requires at most 2.11<5~ 1 quantum queries and this bound 
is in general sharp, see |11|. Hence, for 5 = e/2 we need at most 4.22e _1 quantum queries. 



Obviously, we may reduce the number of quantum queries by choosing a different splitting 
of the two errors for path integration. So, if the first error is, say, ae and the second is (1 — a)e 
for some a G (0, 1), the number of required quantum queries is at most 2.1l£:~ 1 /(l — a), and 
for small a it is roughly at most 2.lle~ 1 . This can be achieved at the expense of increasing 
quantum operations and qubits. 
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We also study the question what is the minimal number of quantum queries for solving 



path integration by an arbitrary quantum algorithm. Similarly as in |L6) , we show that path 
integration is no easier than a specific summation problem. Then using the lower bound of 
I~2|| , we conclude that the minimal number of quantum queries is at least of order e~ 1+a for 



any a G (0, 1), see Theorem 3. This means that the number of quantum queries used by the 
algorithm presented in this paper cannot be significantly improved. 

We stress that the number of qubits is polynomial in e~ l . Furthermore, for the Wiener 
measure the degree is 2 for r = 1, and 1 for r > 2. Hence, if e is relatively large we do not 
need too many qubits to solve path integration on a quantum computer. This is important 
since the number of qubits will be a limiting resource for the foreseeable future. 

From these bounds and from the known complexity bounds in the worst case and ran- 
domized settings, we conclude that 

• Path integration on a quantum computer can be solved roughly e~ l times faster than on 
a classical computer using randomization, and exponentially faster than on a classical 
computer with a worst case assurance. 

• The number of quantum queries is the square root of the number of function values 
needed on a classical computer using randomization. 

We outline the remainder of this paper. In Section 2 we briefly discuss the complexity 
of summation in the worst case and randomized settings, and the quantum summation 
algorithm for computing the arithmetic mean of n numbers, each from the interval [—1,1]. In 
Section 3 we define path integration precisely, while in Section 4 we explain a computational 
approach to path integration. In Section 5 we summarize what is known about the complexity 
of path integration on a classical computer in the worst case and randomized settings. We also 
outline an algorithm of Curbera, ||, which requires exponentially many function values in 
the worst case setting, and which is the basis for the quantum path integration algorithm. In 
Section 6 we discuss path integration on a quantum computer and summarize the advantages 
of the quantum algorithm. In Section 7 we prove that the upper bound on the number of 
quantum queries presented in Section 6 is essentially minimal. In the Appendix we present 
the proof of how many variables must be used to approximate path integrals to within e. 



2 Quantum Summation Algorithm 

Sums occur frequently in scientific computation. For example, when Monte Carlo or Quasi- 
Monte Carlo are used to approximate a <i-dimensional integral, we compute n _1 Ym=i f( x i)i 
where the Xi are ci-dimensional vectors that are chosen randomly (for Monte Carlo) or de- 



terministically (for Quasi-Monte Carlo), see e.g., |13]. As we shall see in Section 5, such 
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algorithms can be also used for approximating path integrals. In fact, for many linear prob- 



lems it is known that linear algorithms enjoy many optimality properties, see e.g., [14], |22], p3| . 
Linear algorithms have the form YH=i a if{ x i) f° r coefficients tij that are sometimes, but not 
always, equal to rT 1 . Let yi = aif(xi)n. Then for all these applications we wish to compute 

n 

S n (y) = n- 1 ^yi. (1) 

8=1 

In this paper we restrict ourselves to the case when \y^\ < 1 for i = 1,2, ... ,n. More 

general conditions on the yi of the form (n~ l Y^i=i \Vi\ P ) — 1 with p G [1, oo] are considered 
in §,0, 

We are interested in applications where n is huge. We wish to approximate S n to within 
e for e G (0, |). The terms yi are not stored or computed in advance. We assume that for a 
given index i we have a subroutine that computes yi. This assumption is typical for scientific 
problems where, as explained above, yi depends on the function value f(xi). 

Before we discuss quantum computation of S n , we briefly mention summation complexity 
results in the worst case and randomized settings on a classical computer, see |14|, [16| . The 
worst case complexity, comp wor (n, e), is defined as the minimal number of operations needed 
to compute an e-approximation to S n for all \yi\ < 1 using deterministic algorithms. The 
randomized complexity, comp ran (n, e), is defined analogously when we permit randomized 
algorithms. It is known that 



and if n ^> e 2 , 



comp wor (ra, e) w n (1 



comp ran (n, e) « e 2 . 



Hence, in the worst case setting we must add essentially all n numbers, whereas in the 
randomized setting it is enough to add only e~ 2 terms and this, of course, can be achieved 
by the Monte Carlo algorithm that chooses e~ 2 samples from the set {yi,y2, ■ ■ ■ ,y n }, each 
with probability n~ l , and computes their arithmetic mean. This speed-up is significant. 

We now turn to what is known about summation on a quantum computer. We wish to 
compute QS n (y,e) which approximates S n (y) to within e with probability at least |. That 
is, QS n (y,e) is a random variable for which the inequality \S n (y) — QS n (y,e)\ < e holds 
with probability at least |. The performance of a quantum algorithm can be summarized 
by the number of quantum queries, quantum operations and qubits which are used, see 
0. ||, ||, [Hj for precise definitions of quantum computation and quantum algorithms. Here, 
we only mention that the quantum algorithm obtains information on the terms y^ by using 
only quantum queries. The number of quantum operations is defined as the total number of 
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bit operations performed by the quantum algorithm, and the number of qubits is defined as 
k if all quantum computations are performed in the Hilbert space of dimension 2 k . 

Since the number of qubits will be a limiting resource for the foreseeable future, it is 
important to seek algorithms which require as few qubits as possible. 

Brassard, H0yer, Mosca and Tapp, see 0, presented a quantum algorithm QS n that 
solves the summation problem. Their algorithm is based on Grover's iterate, see || and 
uses quantum Fourier and Walsh-Hadamard transforms that can be implemented by well 
known quantum gates. Assume that n ^> s^ 1 . Then the algorithm QS n uses of order 

e^ 1 quantum queries, 

e^ 1 log n quantum operations, 

log n qubits. 



More precise bounds are known about the number of quantum queries. In it is shown 



that the quantum algorithm QS n uses at most e 1 2.10 . . . quantum queries and this bound 
is sharp for small e and large en. Due to the lower bound of Nayak and Wu, see Hl2[ , the 



number of quantum queries of any quantum algorithm that solves the summation problem 
must be at least of order e~ x . Hence, the algorithm QS n uses almost the minimal number 
of quantum queries. (In this paper log denotes log 2 .) 

We can run the quantum algorithm QS n several times to increase the probability of 
success. If we want to solve the problem with probability 1 — 6, then we should run QS n 
roughly log 5 _1 times and take the median as our final result. Then the number of queries 
and quantum operations is multiplied by log but the number of qubits stays the same. 

Of course, these quantum results are of interest only if e~ l is significantly less than n. 
Fortunately, this is the case for a number of important problems. Indeed, this paper will 
supply one more such problem, namely, path integration. 

So far we assumed that we summed numbers from the interval [—1,1]. The interval 
[—1,1] is taken only for simplicity. If we have the interval [—M,M] then we can rescale 
the summands to i/i/M, and multiply the computed result by M. This corresponds to the 
previous problem over the interval [—1,1] with e/M. Note, however, that for large M, and 
n > M/e, the quantum cost is of order M/e, which is significantly larger than 1/e. 



3 Definition of Path Integration 

We now define path integrals studied in this paper, see also [j25| . Let X be an infinite 
dimensional separable Banach space equipped with a probability measure /i. We assume 



that /i is a zero mean Gaussian measure, see e.g., [24]. The space X can be embedded in the 
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Hilbert space H = L 2 ([0, 1]) for which the embedding Im : X — > H is a continuous linear 
operator. The inner product of H is denoted by (•, -) H . Then the measure v = /ilm -1 
is also a zero mean Gaussian measure on the Hilbert space H. Let C v be the covariance 
operator of u, i.e., C v : if — > if and 

(C u h 1 ,h 2 ) H = / (h,h 1 ) H (h,h 2 ) H v{dh) Vh 1 ,h 2 eH. 
Jh 

The operator is self adjoint, nonnegative definite and has a finite trace. We can assume 
that there exists an orthonormal system {r/i} from Im(X), (r)i,r)j) H = 5ij, for which 

C u T]i = XiTJi, 

oo 

Ai > A 2 > • • • > and ^ X { < +00. (2) 

i=i 

We illustrate this definition by the important example of the space X = C([0, 1]) of 
continuous functions defined on [0, 1] with the sup norm, \\x\\ = max^^i] The space 

C([0, 1]) is equipped with the classical Wiener measure /i — w. The measure w is a zero 
mean Gaussian measure with covariance function min(t,-u). That is, 



/ x(t)w(dx) = Vie [0,1], 
Jc([o,i]) 

/ x(t) x(u)w(dx) = min(t,-u) Vt, mG[0, 1]. 

•/CfO.l) 



'C([0,1]) 

For the Wiener measure w, we have lm(x) = x and 

'2i - 1 \ . 4 



Vi 



sin f — - — tcx } , \i 



7T 2 (2i - l) 2 ' 



We return to the case of general X and /i. Let F be a class of real-valued /i-integrable 
functions defined on X. The path integration problem is defined as approximating integrals 
of / from F, 

/(/):= / f(x)fjL(dx)= f f(lur l x)u(dx), V/GF. (3) 

If only finitely many eigenvalues A« of C„ are positive, then the measure v is concentrated 
on a finite dimensional subspace of H and path integration reduces to a finite dimensional 
Gaussian integration. To preserve the main feature of the path integration problem, which is 
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integration over an infinite dimensional space, we assume that all eigenvalues A, are positive. 
The element x from H can be written as x = X^i^ 7 ^' with ^ = {x ) r] i ) H . Therefore the 
integrand / in (^[) depends on infinitely many variables U. That is why the path integration 
problem can be viewed as integration of functions having infinitely many variables. 

In this paper we will consider the classes F r of functions whose r — 1 times Frechet 
derivatives exist and are bounded, and whose rth Frechet derivatives satisfy the Lipschitz 
condition. More precisely, for a non-negative integer i, let ||/^|| = sup a . gX 
Here, f^ l '(x) is an i-linear form from X 1 to R, and its norm is defined as = 
su P|k 3 IU<i \f (€> (x)xix 2 ■■■ Xi\. Obviously, ||/ (0) || = ||/|| = sup^ |/(ar)|. 

Let r be a positive integer. For positive numbers Kq, K\ . . . , K r , define ||/|| r -i = 
max <j< r _i The class F r is defined as 

F r = {f : ||/|U<1, \\f ir - 1) (x)-f^ 1 \y)\\ < K r \\1m(x - y)\\ H , Vx,yGX}. 

For r = 1, the class Fi consists of bounded Lipschitz functions. The values of / are 
bounded by Kq, and the Lipschitz constant by K\. For r > 2, the class F r consists of 
bounded smooth functions. All functions from F r are r — 1 times Frechet different iable, their 
ith derivatives are bounded by K{ for % — 0, 1, . . . , r — 1, and the (r — l)st derivatives satisfy 
the Lipschitz condition with the constant K r . 

Note that for any / G F r , the path integral /(/) is well defined since / is continuous 
and bounded. From \f(x)\ < K we have |/(/)| < K . If K is large we can use a different 
estimate on /(/). We have /(/) = f(0) + I(f - /(0)) and |/(/-/(0))| < K 1 J H \\x\\v{dx) < 

( \l/2 / U/2 

K\ \ Y^j=\^3) ■ Hence, |/(/)| < |/(0)| + K\ [YlJLi^j) ■ This estimate can be better 
than the previous one for large Kq. 

As we shall see in the next sections, path integration for the class F r is intractable in 
the worst case setting. We stress that for other classes of functions, path integration can 
be tractable even in the worst case setting. An example is provided for the class of smooth 
integrands occurring in the Feynman-Kac formula, see [HI |T8|| . 



4 Computational Approach to Path Integration 

We want to approximate /(/) to within e for all / G F. The approximate computation 



of /(/) consists of two steps, see [25 . The first is to approximate the infinite dimensional 
integration / by a (^-dimensional integration 1^, where d = d(e, F) is chosen as the minimal 
integer for which the error of this approximation is at most, say, |. The second step is to 
compute an approximation to Id with error at most |. Clearly, we should expect that d(e, F) 
would go to infinity as e goes to zero. 



8 



More precisely we proceed as follows. Let f d : 



be defined for t = [ti, t 2 , . . . , t d ] G 



as 



f d (t) = f (inr^tiTji + t 2 r] 2 + ■■■ + t dVd )) . 



Define 



W) 



1 



1 



(2tt)^ VAiA 2 ---A d 



f d (t)exp(-t 2 l /(2X 1 ) 



(4) 



t 2 d /(2X d )) dt. (5) 



Observe that I d is a finite dimensional Gaussian integral with the eigenvalues Aj as variances. 
Note that the eigenvalues Aj tend to zero. Indeed, since a = A« < +oo and Aj are non- 
increasing then Aj < a/i for all i. Hence, we have decreasing dependence on the successive 
variables £j in (|5|). For continuous /, we have 

/(/) = lim/ d (/ d ). 



As outlined above, we want to choose the minimal d = d(e, F) such that \I(f) — Id(f)\ < § 
V/ G F, and then to compute an |-approximation to a finite-dimensional integral I d (f d ). 
We now find d(e, F r ) for a family of eigenvalues Xj of the covariance operator C„. The family 
includes the eigenvalues of C v for the Wiener measure. 



Theorem 1 Suppose Xj is of order j k with k>l. Then 



cj £ - 2 /( fc - 1 ) < d(e, Fx) < ci e" 
< d(e, F r ) < c r e 



-i/(*-i) 



-2/(fc-l) 

-i/(fc-i) 



/or r > 2, 



•where c r and c r , for r = 1, 2, . . ., are positive numbers independent of e and depending only 
on the global parameters Ki, r, k and the trace Aj- In particular, if Xj = aj~ k with a > 
then 



d{e, Ft) 
d(e, F r ) 



< 1 + 

< 1 + 



(AaKl 



jfe-1 
fc - 1 



i/(fc-i) 



i/(fc-i) 



1 \ 2/(fc-l) 



/or r > 2. 



For the Wiener measure, Xj = 4/(7r 2 (2i - 1) 2 ) ; we have 



7T 2 £ 



1 

tt% 2 



1 

/or r > 2. 
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The proof of this theorem is given in the Appendix. We stress that the upper bounds on 
d(e,F r ) in Theorem 1 depend only on K\ for r = 1, and on K 2 for r > 2, i.e., on the 
Lipschitz constants for / or /', respectively. This means that we can even take all the 
remaining Ki = oo and the upper bounds on d(e, F r ) still hold. On the other hand, the 
lower bounds depend on all of the Ki. The dependence is weak since they only affect the 
multiplicative factors of the power of e' 1 , and the power of e" 1 does not depend on Ki. 



5 Path Integration on a Classical Computer 

In this section we discuss approximation of path integrals on a classical computer. We 
assume the real number model of computation, which is usually used for the analysis of 
scientific computing problems, see |21| for the rationale. We assume, in particular, that 



we can perform arithmetic operations (addition, subtraction, multiplication, division), and 
comparisons of real numbers. We assume that these operations are performed exactly and 
each costs unity. To approximate path integrals we must have information concerning the 
integrands / G F. This information may be supplied by function values f{x,j) for some 
x\, x 2 , . . • , x n , where n = n(e, F) will be chosen depending on the error demand e and the 
class F. As outlined in the previous section, we will need to know f(x) for x belonging 
to a finite dimensional subspace X d = span(Im _1 ^i, Im _1 ?72, . . . , Im -1 ^) with d = d(e,F). 
We therefore assume that we can compute values of f(x) for x G X d and the cost of one 
such evaluation is c d . Usually ^> 1. Furthermore, we will sometimes assume that the 
cost c d depends linearly on d, i.e., c d = cd; however, this assumption is not essential to 
the analysis. For a more complete discussion of the real number model of computation with 
function values, see |L5], |2"2"| . 

Let A(f) be any algorithm for approximation of path integrals. The algorithm A uses 
a finite number rii of function values at points Xi and a finite number n 2 of arithmetic 
operations and comparisons to compute A(f). The cost of computing A(f) is c d rii + n 2 . 
In the worst case setting, the error and cost of A are defined by its worst performance over 
the class F . In the randomized setting, the algorithm A may use randomly chosen samples 
Xi, and its error and cost are defined by the expected error with respect to the distribution 
generating the random samples for a worst / from F. By the worst case or randomized 
complexity, we mean the minimal cost that is needed to compute an ^-approximation for all 
/ G F, see |22j for precise definitions. 



We now briefly discuss the worst case and randomized complexities of path integration 
for the classes F r . We begin with the worst case setting. We first state the result of 



Bakhvalov, see e.g., |14| , |22| , which states that the worst case complexity of multivariate 
integration over the unit cube [0, l] d for r-times differentiable functions is of order c d e~ d / r . 
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For path integration d is an increasing function of e~ x , and as shown in Theorem 1, it goes 
to infinity polynomially in e~ l as e goes zero. This suggest that the worst case complexity, 
comp wor (e, F r ), of path integration in the class F r is exponential^] in e~ l . A formal proof may 
be found in |25| for any r, and more precise complexity bounds in || for r — 1. Thus path 
integration is intractable for the class F r in the worst case setting. This means that the cost 
of any algorithm for solving this problem must be exponential. Yet, as we shall see, such 
algorithms will be useful for quantum computation. We now sketch such an algorithm. 

We first consider the case r = 1 and then show that an easy modification of the same 
algorithm can be also used for r > 2. We assume^ that Xj = B(j~ fe ). From Theorem 1 we 
know that it is enough to compute an |-approximation to the integral Id(f) with d = d(e, Fi) 
given in Theorem 1. We have 

W) = I fd(t)v d (dt), (6) 

where v& is a Gaussian measure on M d with mean zero and with the diagonal covariance 
matrix diag(Ai, A 2 , . . . , \d). 

This problem has been studied in ||. Based on that paper we describe an algorithm S n 
with worst case error at most |. We opt here for simplicity of the presentation of S n at a 
slight expense of its cost. Let n = m d for the minimal odd integer m for which 

m > ^ — — . (7) 

For x > 0, let ip(x) = J q x exp(— t 2 /2)dt be the probability integral, and let ip~ l be its 

inverse. We note that it is easy to compute numerically for any t e R. As in Lemma 1 

of m, for i = 1, 2, . . . , d define the points tij : 



U,o = -oo, 

t itj = (3A i ) 1 / 2 ^- 1 ((i-l/2)/m), j = 1,2, 

tim+l = °°- 



m, 



Then take t*j = t^j if \t^j\ < and t* • = t it j + i otherwise. For the integer vector 

j = j 2 , . . . ,jd], with ji = 1,2, ... ,m, define the sample points 



x 



: We follow a convention of complexity theory that if the complexity growth is faster than polynomial 
then we say it is exponential. 

2 The theta notation means that there exist positive numbers c\ and C2 such that c\j~ k < Xj < C2]~ k for 
allj = 1,2,.... 
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Then the algorithm takes the simple form 



S n (fd) 



= n 



8 



3 



Curbera proved in [[3] that the worst case error of S n is at most 



d 



) 



1/2 



e WOI (S n ) <2K 1 [n 



m 



i 



£ 

< - 

~ 2 



where the last inequality holds due to the choice of m. 
The cost of S n is (c<f + l)n where 




=r e -a/(*-i) 



n = n(e 



(9) 



where a = 4Ki (tt Yli=i ^i) + 2e, and cT is from Theorem 1. 

We now consider the case r > 2. As shown in Theorem 1, we can now restrict ourselves 
to the integrals Id(f) for d = d(e,F r ). We stress that d(e,F r ) is much less than d(e,Fi) 
for small e. Observe that all functions / from F r also belong to Fi since they satisfy the 
Lipschitz condition with the constant K\. Hence we can use the algorithm S n with the 
important difference that now d = d(e,F r ). Hence, we compute an |-approximation by the 
algorithm S n with cost (c<i + l)n, where 



with C2 from Theorem 1. 

We now justify why it is enough to apply the algorithm S n for the class F r for any r > 2. 
The reason is that for path integration the smoothness parameter r is not as important 
as for finite dimensional integration. Indeed, since the exponent of e -1 for the worst case 
complexity of path integration is unbounded (as e — > 0) for any fixed r, it does not help 
much to divide by r. 

As we shall see in the next section, for quantum computation the logarithm of the worst 
case complexity is important and r can only effect a multiplicative factor. The most im- 
portant property is how fast d(e,F r ) goes to infinity. As we know from Theorem 1, the 
influence of r is significant here since we have different formulas for d(e,F r ) for r = 1 and 
r > 2. However, for r > 2, the use of more efficient algorithms than S n can only improve 
the multiplicative factor of the logarithm of the worst case complexity. 




(10) 
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We turn to the randomized setting for Aj = 0(j~ fc ). The randomized complexity, 
comp ran (£:, F r ), can be easily obtained by applying results of Bakhvalov for finite dimen- 
sional integration. The analysis in |25| yields 



comp ran (£, F r ) = 6 ((c d + 1) e~ 2{1+o{1)) ) as e -> 0, (11) 

where d = d(e, F r ) and the factors in the notation depend at most quadratically on Kq 
and K\. Hence, in the randomized setting we have roughly quadratic dependence in e~ l on 
the number of function values. If = c d then 

f /l\ i T 2 f i + 2+0 « x 
comp ran (e, F r ) = 6 c - j (12) 

where 7(1) = 1 and 7(7*) = for r > 2. 

Hence, the randomized complexity of path integration depends polynomially on e^ 1 , 
and therefore the path integration problem is tractable in the randomized setting. In fact, 
the upper bound can be achieved by the Monte Carlo algorithm with randomized error at 
most I and with the cost proportional to KqE~ 2 randomized evaluations of a function of 
d = d(e, F r ) variables, where d(e, F r ) is given by Theorem 1. 

Note, however, that if k goes to one then the degree of e~ l in the randomized complexity 
goes to infinity. The reason is that in this case we have to compute function values of very 
many variables. On the other hand, for the Wiener measure we have k = 2, and the degree 
of e' 1 is roughly 3 + j(r). 



6 Path Integration on a Quantum Computer 

We now analyze path integration on a quantum computer. The idea behind solving path 
integration on a quantum computer is quite simple. (However, the analysis is not so simple.) 
We will apply analogous techniques for other problems in future papers. 

Start with an algorithm that computes an ^-approximation to path integration in the 
worst case setting and that requires summation of the form of ([!]). We run this algorithm 
on a quantum computer using the quantum summation algorithm of Section 2. Obviously, 
n is now a function of e~ l . For path integration for the class F r we know that n is an 
exponential function of e~ l and is bounded by @ for r = 1, and by ( |10D for r > 2. However, 
the exponential dependence on e~ x is now not so essential since the cost of the quantum 
summation algorithm QS n depends only logarithmically on n. Since log n is a polynomial 
in e~ l we conclude that path integration on a quantum computer can be solved at cost 
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polynomial in e~ l . That is, intractability of path integration in the worst case setting is 
broken on a quantum computer by the use of the quantum summation algorithm. 

For other intractable problems in the worst case setting for which the worst case com- 
plexity can be achieved by summation of n numbers, intractability will be broken as long 
as n is a single exponential function of i.e., n (e- 1 ) < 2 p ( £ ) with p being a polynomial. 
Then the quantum cost will be polynomial in and the problem will be tractable on a 
quantum computer. This idea will not work if n(e~ l ) is a double exponential function (or 
worse) of e~ l since then the logarithm of n^ 1 ) will be still an exponential function of e^ 1 . 

We now provide details of this idea for path integration for the class F r with eigenvalues 
Xj = Q(j~ h ). We take the algorithm S n defined by (^) with n given by (^) for r = 1 and 
by fllOD for r > 2. The algorithm S n already has the form (H) required by the summation 
algorithm. However, the summands f(x-j) are not necessarily in the interval [—1,1]. The 
function / belongs to F r and therefore its values are bounded by Kq. Hence, it is enough to 
scale the problem by running the quantum summation algorithm for yj = f(xj)/Ko, replace 
e by e/ Kq, and multiply the computed result by K . The cost of an algorithm on a quantum 
computer using m qubits is defined as on a classical computer with the cost of a quantum 
query taken asc^ + m since /(xj)'s are computed and m qubits are processed by a quantum 
query. Using the results of quantum summation from Section 2 applied for S n we obtain the 
following theorem. 

Theorem 2 Consider path integration for the class F r with the eigenvalues Xj = 0(j~ fc ). 
Using the quantum summation algorithm QS n to compute an e / K - approximation to S n , we 
compute an e- approximation for path integrals with probability at least | and of order 

e^ 1 quantum queries, 

£ -{(k+j{r))/(k-i) j Q g g -i quantum operations, 

e -((l+7(r})/(M) bg e -l qubitS} 

where 7(1) = 1 and j(r) = for r > 2. If = cd, then the cost of this algorithm is of 
order 

l + 7 (r) 1 

i\ + V _ L1 1 

c + log - 



ej \ e t 

For the Wiener measure the results are more precise. The algorithm requires (neglecting 
ceilings for simplicity) at most 

2KqE V quantum queries, 
a p {e,r r )[og quantum operations, 

d up (e,F r ) log 16K ° Kl quhl ts. 
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If Crf = c d, then the cost of this algorithm is at most 



2K d up (e,F r 



£ 



c + 2 log 



16 K K 



£ 



) 



Here d up (£, F r ) is an upper bound on d(£, F r ) given by 



ir 2 £ 



) 



l+ 7 (r) 




^(e,^) 



with (5{r) = 2Ki for r = 1, and (3{r) = K 2 for r > 2. 

If we want to increase the probability of computing an ^-approximation to path integration 
then, as explained in Section 2, we can run the quantum algorithm for QS n roughly log 5~ l 
times and take the median as the final result. Then the probability of success is at least 
1—5. Obviously the cost is then multiplied by log <5 _1 but the number of qubits stays the 
same. 

We compare cost(QS n ), the cost of the quantum algorithm, with the worst case complex- 
ity of path integration. The essence of Theorem 2 is that cost (QS n ) depends polynomially on 
£~ l . Since the worst case complexity is exponential in the use of quantum summation 
breaks intractability of the worst case setting. Note that we have exponential speed-up, i.e., 
comp wor (£:, F r )/cost(QS n ) is exponential in e~ x . 

We now compare cost(QS n ) with the randomized complexity of path integration. As dis- 
cussed in Section 5, path integration is tractable in the randomized setting and its randomized 
complexity is characterized by (|TT|) and fll2"|). Comparing the formulas for the randomized 
complexity with cost(Q5' n ) we see that the ratio of the number of quantum queries used 
by the quantum algorithms to the number of function values used by the best randomized 
algorithm is roughly £ . If we compare cost(QS n ) to the randomized complexity we see 
that the speed-up is roughly of order e -1 . That is, we solve path integration on a quantum 
computer roughly e _1 times cheaper than on a classical computer using randomization. We 
summarize our results in the following corollary. 

Corollary 1 Consider path integration for the class F r with Xj = ©(j~ fc ). Then 

• Path integration on a quantum computer is tractable. 

• Path integration on a quantum computer can be solved roughly e~ Y times faster than on 
a classical computer using randomization, and exponentially faster than on a classical 
computer with a worst case assurance. 
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• The number of quantum queries is the square root of the number of function values 
needed on a classical computer using randomization. 

• The number of qubits is polynomial in e~ x . Furthermore, for the Wiener measure the 
degree is 2 for r = 1, and 1 for r > 2. 



7 Lower Bounds on the Number of Quantum Queries 

We now study lower bounds on the minimal number qq(e, F r ) of quantum queries needed 
to compute an e- approximation with probability 3/4 for path integration for the class F r . 
From Theorem 2 we know that qq(e, F r ) is at most of order e~ x . We show that this bound 
cannot be significantly improved. 

Theorem 3 Consider path integration for the class F r with all positive eigenvalues \ . Then 

lim e l ~ a qq{e, F r ) = oo VaG(0,l). 

Proof: The proof consists of two steps. The first one is to reduce path integration to a 
finite dimensional Gaussian integration which is no harder than the original problem. The 



second step is essentially the same as in Novak's papers, see [[L4|, |16| , and reduces the finite 
dimensional Gaussian integration problem to summation for which the lower bound of Nayak 
and Wu, see [0, applies. 

In the first step of the proof, for a given a G (0, 1) we take an integer d > r(l — a) /a. 
(Hence, d is large for small a.) The path integration problem for the class F r is no harder if 
we assume some additional properties of functions / from F r . We have, see (|3]), 

/(/) = [ film^x) u{dx) 

H 



where x = YliLi UVi- 

Let us now assume that / depends only on the first d components t\, t 2 , . . . , td, and call 
this class F r d . Obviously, F r)d is a subclass of F r and therefore path integration for F r>d is 
no harder than for the class F r . For the class F rjd , the path integration problem reduces to 
a finite dimensional Gaussian integration problem. That is, for / e F r ^ we have 



W) = / fd(t)u d (dt), 
where v d is the Gaussian measure given by ([]), and f d is given by (||). 
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The functions f d from F r>d are r — 1 times continuously differentiable and 



H/rfll = sup |/ d (t)| < K . 



Furthermore, their r — 1 partial derivatives satisfy the Lipschitz condition. More precisely, 
there exists a positive number (3\ depending only on d, r and K Q , Ki, . . . , K r such that 

|i>7d(0--D7d(y)| < AP-j/IL vt, y eR d , 



where D l runs through the set of all partial derivatives of order r — 1. 

d 



This shows that the class F ryd is closely related to the class 1,1 studied by Novak, see 



0, 

Fr 1 ' 1 = {f:[0,l] d ^R\feC r -\[0,l] d ), \D<f(t)-D*f(y)\ < \\t - yU Vt, y }. 

Obviously, the different Lipschitz constants: /?i in our case and 1 for the class F d ~ x ' do 
not play a major role since they do not change the order of error bounds. One difference 
between the two classes is that the common domain of functions from F r ^ is ~R d , whereas for 
the class -F^ -1 ' 1 the common domain is [0, l} d . A second difference is that we have Gaussian 
integration whereas Novak considered uniform integration, J, Q1 , d f(t) dt. As we shall see 
below these two differences are not really essential. 

In the second step of the proof, we use Novak's proof technique. From |T3, |TB| we know 
that for any positive E\ there are functions fi, f 2 , ■ ■ ■ f n , with n = Q(e l d ^ r ), from the class 
Fj -1 ' 1 such that they take non-negative values, have disjoint supports in [0, l] d and 



L 1]d fi(t)dt = e\ +d/r i = l,2 



Uo,i] 



n, 



• Y^=l a ifi e F d ' V N < L 

We use the same functions fi for our Gaussian integration for the class F r ^- Since the 
support of fi is in [0, l] d we can extend fi by zero to M. d . The extended functions fi have 
exactly the same smoothness as required for the class F r ^ d , and there exists a positive /3 2 
depending only on d, r and K , K lf . . . , K r but independent of i such that j3 2 fi £ F r>d . Note 
that 



fi(t)u d (dt) = / p d (t) fiit) dt, 

J[0,l] d 

where 
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Since all Xj are positive, the function has positive minimum and maximum over [0, l] d . 
That is, there are positive /3 3 and /5 4 depending on d and Ai, A 2 , . . . , A^ such that 



ft < Pdit) </? 4 Vt,- G [0,1]. 
Therefore for the functions gi — /3 2 fi G F r( j we have 



(3 2 {3 3 e l+d / r < I( 9i ) = [ g t (t)u d (dt) < (3 2 f3 A e\ +d/r ', (13) 



and 



Since 



^ cKiffj G F r , d V|a»| < 1. 



n 



\i=l / i=l 

we reduce our problem to summation of n — 0(6^ r ) terms for arbitrary \an\ < 1. Let 



/3 2 f3 i e 1 



l+d/r 



Oti. 



Observe that by varying oti from [—1, 1], the y^ can take any value from /3 3 //3 4 [— 1, 1] due to 
the left hand side of flT3[). 

We need to compute an ^-approximation to I(X^Li ^ift)- This is equivalent to computing 
an ^-approximation to 

1 n 



n 
1=1 



with e 2 = e/(n(3 2 (3 4: e 1 1 +d ^ r ) = Q(s/si). 

The summation problem n _1 Y27=i Vi f° r our Vi ^ s n °t eas i er than the summation problem 
n ~ XS LA=\Vi f° r — fls / 'Pi- We can now apply the lower bound of Nayak and Wu, 

see Jl2| , that states that the minimal number of quantum queries needed to compute an e 2 - 
approximation with probability 3/4 for the summation problem n~ l Y^i=i y% with \yj\ < (3%/ (3$ 
is bounded from below by 



C min \ n, = © (nrinfo d/r , e x /e) ) , 



with some absolute positive number C . 
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Finally, we take e\ such that e\ d = £i/e, i.e., S\ = e r ^ d+r \ and conclude that the 
minimal number qq(e, F r ) of quantum queries is at least of order £~( 1 ~ r /( d + r )). Since d > 
r(l — a) /a implies that a > r/(d + r), we have 

e l - a qq(e, F r ) = Q ( £ -(<x-r/(d+r))j _^ ^ as e ^ . 
This completes the proof. ■ 



8 Appendix 

We prove Theorem 1. We begin with r = 1. It is shown in [^5| that d = d(e, Fi) < d* where 
d* is an integer satisfying 



A, < e 2 /(2K 1 



i=d*+i 



For Aj = Q(i ) with k > 1, we get 



d* = 6 ((K^e) 2 '^) as e -> 0. 



For A, = az , we get 



a \ 1/(W) /2^^ /(M 



k — 1 



For the Wiener measure we have 



1 2K t y 1 
+ 2 



7T 2 £ 



This establishes upper bounds on d(e, F\). 

To get a lower bound, take the function g : R — ► R defined by g(x) = ci\x\/(l + \x\) with 
ci = mxn(K , Kx). We have sup,j, eR \g(x)\ — C\ < K , and g satisfies the Lipschitz condition 
with the constant C\ < K 1 . For x G H, define x d = ^j=d+i ( x iVj) H^i an< ^ 

f{\m l x) = g{\\x d \\). 

Then f & Fi and / d = 0. We have 1(f) - I(f d ) = 1(f) and 



g(\\x d \\) v(dx) > d 



|a;d||z/((ix) 



XJ <i 
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Since (j2% d+1 a f) ' > d 1/2 \ a i\ for an Y a i e R > we § et 

2d 

'(f) > ^ £ K 



where 



d „ / 2d \ 

Intj = JJ(27rA,)" 1/2 / exp - ^ x 2 /(2A;) J da;. 

i=l JYjltd+i x 1< 1 \ i=d+l J 

There exist two positive numbers a% and a 2 such that a.\i~ k < Aj < a 2 i~ k for all z. By 
changing variables = ^/(A;) 1 / 2 and noting that a 2 d~ k > A^ > A« > A2d > ai(2d)~ k we 
conclude that 

Int, = ^^J V ,|expf-^t 2 /2Ux, 



(2vr) d / 2 

Here, t = [ti, £2, • • • , td] and ||£|| = (X^=i t 2 ) 1 ^ 2 . The first integral is just y/2/n. We now 
show that the second integral goes to zero with d. Indeed, let be for the Gaussian measure 
on W 1 with zero mean and the identity covariance operator, and let a = (d k / 'ol 2 ) x I 2 . Then 
the second integral is J Rd \ B \ti\v d (dt), where B a denotes the ball of radius a, and is not 
greater than 

t 2 Mdt)) {I - u d {B a )f 2 . 



The integral with the integrand t\ is just one, and using Lemma 2.9.2 from |22[] p. 469 we 
conclude that 

l-P d {B a ) < 5exp(-a 2 /(2c/)). 

Since a 2 / (2d) = d k ~ l a 2 /2 and k > 1 then this ratio goes to infinity, and 1 — z/ d (i? d ) goes to 
zero. This means that Intj is at least of order d~ k ^ 2 , and /(/) is at least of order d - ^" 1 ^ 2 . 
Hence to guarantee that /(/) = /(/) — I(fd) < § we must take d of order £?~ 2 /( fe_1 ) which 
completes the proof for the case r = 1. 

Assume now that r > 2. We first establish an upper bound on d(e,F r ). For x = 
YlJLi'tjVj £ H, define t = [ti, t 2 , ■ ■ ■] G M°° and t d = [f 1; t 2 , ■ ■ ■ , td]. Then we can identify 
fit) with /(Im _1 x) and f(t d ) with fd(t). By Taylor's theorem we have 

/(f) = f(t d ) + f'(t d )(t - t d ) + / (/' {t d + u(t - t d )) - f{t d )) (t - t d ) du. 
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Note that t — t d = [0, . . . , 0, t d+ i, t d+2 , ■ ■ .] and since f'(t ) is a linear form we have 



f'(t d )(t-t d ) = 

where aj = aj(ti, t 2 , • • • , td)- The mean element of v is zero, which implies that I(f'(t d )(t - 
t d )) = 0. Hence 

1(f) - I{h) = I ( jf 1 (/' {t d + «(* " t d )) - f(t d )) (t - t d ) du^ . 
For r > 2, /' satisfies the Lipschitz condition and we get 

W)-W)\ < K 2 l(\\t-t d \\ 2 J\dv)j 



Vj=d+1 



#9 



For Aj = 6(j fc ) we obtain 

= o(d-<*- 1 )) 

and for d = 0(e-V(*-i)) we guarantee that |J(/) - J(/ d )| < | for all / G F r . 



For A,- = aj k , we have 



j=d+i 

In this case it is enough to take 

d = 



°° POD 

^ Aj < a / u" k du 



2(fc-l) d*- 1 " 



fc - 1 



V(fc-i) / 1 \V(fe-i) 



For the Wiener measure, 

j:x,<- 2 j d («-i/2)- 2 * = ^ 

and 



1/2)' 



^ ll 

7T% 2 
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This establishes upper bounds on d(e, F r ). 

To get a lower bound, consider the function g(x) = c r x 2 / (1 + x 2 ) = c r (l — 1/(1 + x 2 )) for 
x G R, where c r is a positive number chosen such that maxo<i< r _i sup^g \g^(x)\/Ki < 1, 
and such that g^" 1 ' satisfies the Lipschitz condition with the constant K r . It is easy to see 
that such a positive number c r exists. Indeed, the j'th derivatives of 1/(1 + x 2 ) can be written 
as the ratio of two polynomials pj(x)/(l + x 2 y +1 with the degree of pj being at most j, and 
therefore all derivatives go to zero as \x\ goes to infinity. 

As for the case r = 1, we take Xd = Y^f= d +\ ( x > ^i) h Vj: an d /(Im _1 x) = g(||xd||)- Then 
/ G F r , f d = and /(/) — I(fd) = 1(f)- Similarly as for r = 1 we have 

Hf) > ^ / \\t\\ 2 u d (dt), 

J||t||<a 

where a = (d k / a 2 ) 1 ^ 2 . We now show that the last integral tends to d. Indeed, it can be 
written as 

\t\\ 2 u d (dt) - l \\tfv d (dt). 

J \\t\\>a 

The first integral is obviously d, and we show that the integral over the outside of the 
ball tends to zero. For large d, the norm of t is also large, and we can estimate ||t|| 2 < 
exp(Q||t|| 2 /2) for c d = 2a 2 \n(d k /a 2 )/d k . Then 

t\\ 2 u d (dt) < (l-c d ) d/2 [ u diC (dt), 

t\\>a J\\t\\>a 

where Vd )C is a Gaussian measure on IR d with mean zero and covariance operator (1 — c d )~ x I. 



Again using Lemma 2.9.2 from [22] we obtain 



\tfv d (dt) < 5(1 -c d ) d/2 exp (-a 2 (1 -c d )/d). 



t\\>a 



Since k > 1, the quantity (1 — c d ) tends to 1, and since a 2 /d = Q(d k ~ 1 ) tends to infinity, 
the integral goes to zero as claimed. 

Hence, /(/) is at least of order d^^ 1 ^ and d must be at least of order e -1 /^" 1 ) to 
guarantee /(/) < |. This completes the proof for r > 2. ■. 
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